Magnetism in graphene nano-islands 
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We study the magnetic properties of nanometer-sized graphene structures with triangular and 
hexagonal shapes terminated by zigzag edges. We discuss how the shape of the island, the imbalance 
in the number of atoms belonging to the two graphene sublattices, the existence of zero-energy states, 
and the total and local magnetic moment are intimately related. We consider electronic interactions 
both in a mean-field approximation of the one-orbital Hubbard model and with density functional 
calculations. Both descriptions yield values for the ground state total spin, S, consistent with Lieb's 
theorem for bipartite lattices. Triangles have a finite S for all sizes whereas hexagons have 5 = 
and develop local moments above a critical size of « 1.5 nm. 

PACS numbers: 



The study of graphene-based field effect devices lias 
opened a new research venue in nanoelectronicaii^ii^ii^i^. 
Graphene is a truly two-dimensional zero-gap semicon- 
ductor with peculiar transport and magnetotransport 
properties, including room temperature Quantum Hall 
effeclj^, that makes it very different from conventional 
semiconductors and metals^. Progress in the fabri- 
cation of graphene-based lower dimensional structures 
have been reported both in the form of one-dimensional 
ribbons^i^ and zero-dimensional dots^ii^ii^. Interestingly, 
the electronic properties of graphene change in a non- 
trivial manner when going to lower dimensions. Ribbons, 
for instance, can be either semiconducting with a size de- 
pendent gap or metallio^i^. 

The electronic structure of graphene-based nanostruc- 
tures is expected to be different from bulk graphene 
because of surface, or, more properly, edge effectsi^. 
This is particularly true in the case of structures with 
ziz-zag edges which present magnetic propertiesi^'i^'i^. 
Whereas bulk graphene is a diamagnetic semimetal, sim- 
ple tight-binding models predict that one-dimensional 
ribbons with zigzag edges have two flat bands at the 
Fermi energy^^'^^'— li^ii^, i.e., are paramagnetic metals. 
Spin polarized density functional theory (DFT)ii and 
mean field^ calculations confirm that these bands are 
prone to magnetic instabilities. 

The fabrication of graphene nanostructures using top- 
bottom techniques does not permit creating atomically 
defined edges to dateiS. In contrast, bottom- up pro- 
cessing of graphene nano-islands is very promisingi^. 
Hexagonal shape nano-islands with well-defined zigzag 
edges atop the 0001 surface of Ru have already been 
achieved^. This experimental progress motivates our 
study of the electronic structure of graphene nanostruc- 
tures with various shapes. Graphene quantum dots also 
hold the promise of extremely long spin relaxation and 
decoherence time because of the very small spin-orbit and 
hyperfine coupling in carbor^ii. 

We have found that, remarkably, both the DFT calcu- 
lations and the mean field approximation of the single- 



band Hubbard model with first-neighbors hopping yield 
very similar results in all cases considered. Our mean 
field results are in agreement with the predictions of 
Lieb's theoreni^i regarding the total spin S of the ex- 
act ground state of the Hubbard model in bipartite lat- 
tices. The honeycomb lattice of graphene is formed by 
two triangular interpenetrating sublattices, A and B. 
Triangular nanostructures have more atoms in one sub- 
lattice, Np^ ^ Nb; our mean field calculations consis- 
tently predict that the total spin of the ground state is 
25" = Na — Nb and that is mainly localized on the edges. 
This could have been anticipated from Hund's rule and 
the appearance, in the non-interacting model, of Na—Nb 
degenerate states with strictly zero energy. Hexagonal 
nanostructures, for which Na — Nb, result in 5 = 
ground states even when interactions are turned on. A 
value of = does not preclude, however, an interesting 
magnetic behavior. In fact, we predict a quantum phase 
transition for hexagons: Whereas small ones are param- 
agnetic, large ones are compensated ferrimagnets, both 
with 5* 0. 

The shape and the single-particle spectrum.- The dif- 
ferent atomic structure of triangular and hexagonal 
graphene nanostructures can be appreciated in Figs, [^a) 
and (b). Zigzag edges are formed by atoms that belong 
to the same sublattice, A or B. In the case of the tri- 
angular systems the three edges belong to the same sub- 
lattice, hereafter A, whereas in the hexagon three edges 
are A-type and the other three are B-type. The edge im- 
balance in triangular nanostructures results in a global 
imbalance, so that the total number of atoms in the sub- 
lattice A and B is not the same. In what follows we 
characterize the size of both triangular and hexagonal 
nanostructures by the integer number N of edge atoms 
of the same sublattice along one edge of the island (see 
Fig. [J). 

The structural differences between hexagonal and tri- 
angular nanostructures reflect in their electronic prop- 
erties. At the simplest level, we describe them in the 
one-orbital tight-binding approximatio n^^i^'^d^d^d^ . The 
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FIG. 1: (Color online), (a) and (b) Atomic structure of 
the triangular and hexagonal graphene islands. (c),(d) Single 
particle spectra for the N = 8 triangle (left) and hexagon 
(right) (e),(f) Sublattice resolved edge content (eq. [1} and 
sublattice polarization (eq. [2| 



model Hamiltonian Hq is totally defined by the positions 
of the atoms and the first-neighbour hopping parameter 
t, which we take equal to 2.5 eV. We set the on-site en- 
ergies for all the carbon atoms equal to zero. We assume 
that the edge atoms are pasivated, so that there are no 
dangling bonds. In Figs. [Ijc) and (d) we show the (low 
energy) spectra corresponding to a triangle (left) and a 
hexagon (right) with edge size = 8. If the system is 
charge neutral the relevant electronic states, correspond- 
ing to the highest occupied and lowest unoccupied molec- 
ular orbitals (HOMOs and LUMOs) are around E = 0. 
The most striking difference between the spectra of the 
triangle and the hexagon is the existence of a cluster of 
zero-energy states in the case of the triangle. A suffi- 
cient condition to have Nz states with strict zero energy 
in graphene structures is to have a sub-lattice imbalance 
Nz = Na — Nb- In the case of graphene islands with 
triangle shape (see Fig. [T]), the sublattice imbalance sat- 
isfies Na- Nb = N -I. 

In order to quantify the edge/bulk character of the sin- 
gle particle eigenstates 0„ (/) , we define their sub-lattice 
resolved edge content: 



I^rj^edge 



(1) 



where rj — A, B and / runs over the N^j atoms. We also 
define the sublattice polarization 

>cH-Ei'^»wi'-Ei'^"Wi' (2) 

leA leB 

Both in the case of the triangle [Fig. [T] (e)] and the 
hexagon [Fig. [1] (f)] states have a predominant edge 



character close to the Dirac point {E = 0), but, again, 
there are some differences. In the triangle, the zero en- 
ergy states have a full sublattice polarization £ = 1 and 
their edge content can reach almost 1. In the case of the 
hexagon there is a perfect AB symmetry (C — 0.5) and 
the edge content does not go above 0.8. 

Electron- electron interactions.- The manifold of 2Nz 
zero energy states, including the spin, of the triangle is 
half-filled. Electronic repulsions determine which of the 
2^^ spin configurations has the the lowest energy. If 
Hund's rule operates in this system, the ground state of 
triangular graphene nanostructures (or any other sub- 
lattice unbalanced graphene systems, for that matter) 
should have a maximal magnetic moment 2S = Nz- In 
contrast, the single-particle spectra of hexagons features 
some dispersion, which acts against interaction induced 
spin polarization. To put this on a quantitative basis, 
we have calculated the electronic structure using both a 
mean field decoupling of the one-orbital Hubbard model 
and DFT calculations in a generalized gradient approx- 
imation (GGA) as implemented in the GAUSSIAN03 
code^^, using an optimized minimal basis seli^^. 
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FIG. 2: (Color online). Left column: self-consistent energy 
spectra for a graphene triangular island with N = 8 (fig l.a). 
Closed (empty) symbols correspond to full (empty) single par- 
ticle states. Right Column: local magnetization close to one 
of the corners of the triangle. Uper row: DFT calculations. 
Lower Row: mean field calculations with the Hubbard model. 
Magnetization arrows are plotted horizontally for the sake of 
clarity 

In the mean field approximation for the Hubbard 
model we solve iteratively the Hamiltonian 



H = Hq + U ^ni^{nii) + nii{ni^), 



(3) 



where Hq is the single particle Hamiltonian described 
above and (n/o-) is the statistical expectation value of 
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the spin-resolved density on atom /, obtained using the 
eigenvectors of H. This mean field decoupling can de- 
scribe spontaneous symmetry breaking along a chosen 
axis. The results shown here were obtained fixing 7V| 
and Ni, with N-^ + equal to the number of carbon 
atoms in the structure. This permits to compare with 
DFT calculations where one typically fixes iVf — Ni. 
A self-consistent solution of H is characterized by the 
spin density: m/ = !LLrz!LLi^ and a single particle spec- 
trum e„^CT. The total spin S = mj obviously satisfies 

^ - 2 ■ 

Uncompensated lattices: Triangles.- In Fig. [2]we show 
the spectrum and the spin density for a = 8 triangle. 
Upper panels correspond to DFT results with hydrogen 
passivation of the edge atoms. The results in the lower 
panel correspond to the mean field results for the Hub- 
bard model. In both cases we have verified that the so- 
lutions with iV-f — TVj^ = Nz = 7 have lower ground state 
energy than solutions with different value of 25'. The 
typical energy differences are above 0.5 eV. We choose 
the value of U such that the HOMO-LUMO gap in the 
majority spectrum is the same. In the case shown in Fig. 
[2]this corresponds to U = 3.85 eV. Notice that the mean 
field and DFT spectrum have very similar structure in 
the neighbourhood of Ep. Interactions open a spin gap 
in the single-particle zero-energy manifold, resulting in 
maximal spin polarization of those states. The magne- 
tization density of both calculations is also very similar: 
The A atoms on the edge are copolarized positively (right 
arrows) and their B neighbour atoms are copolarized neg- 
atively. The net total spin is mostly sitting on the edge 
and the local magnetization goes to zero in the center of 
the island. Using the same procedure as above to fix [/, 
we find that its value decreases as the size of the islands 
increase. The values of U so obtained are always below 
the critical value U ~ 2.2t ~ 5.5eU above which infinite 
graphene would become antiferromagneti o^^'^'^ . 

These results indicate that the Hubbard model cap- 
tures the low-energy physics of graphene triangular nano- 
islands. One can conclude that next-to- nearest neighbor 
hopping, long-range Coulomb interactions, and correla- 
tions, as included in the DFT calculations, have a minor 
effect on the low energy sector. Importantly, the basic 
features of the mean field solutions, like the structure of 
the spectrum, the total spin of the ground state and the 
magnetization density, are very robust with respect to 
the value of U. We have found very similar results for 
triangles with N between 5 and 30. The solution that 
minimizes the ground state energy always satisfies 

2S - Ni = Na - Nb = Nz = N - 1. (4) 

Our mean field Hubbard model and DFT results are 
in agreement with the Lieb theorem that states that the 
spin S of the ground state of a Hubbard model in a bi- 
partite lattice satisfies the relation 25* — Na ~ Nb^^- 
If the Hubbard model with first-neighbors hopping and 
constant U can be applied to graphene-based structures 
of arbitrary shape, the theorem permits to predict the 



spin of the ground state by simple counting of the sub- 
lattice imbalance. The fact that the number of strict 
zero energy states Nz equals to Na ~ Nb provides a sim- 
ple picture of how the magnetization comes about: Spin 
polarization results from Hund's rule and the absence of 
kinetic energy penalty in sublattice unbalanced graphene 
structures. 

Compensated lattices: Hexagons.- In the case of bal- 
anced structures Lieb's theorem predicts that they have 
5 = 0. This is compatible with a locally unpolarized 
state, but also with locally polarized solutions with an- 
tiferromagnetic correlations. In these cases, calculations 
are necessary to obtain the local magnetization density. 
In the case of hexagons there is a competition between 
the dispersion of the single-particle spectra and interac- 
tions. Dispersion occurs because of the hybridization of 
states that otherwise would lie in a single sublattice close 
to the edge. These states overlap in the inner region and 
close to the vertices and hybridize through hopping in Hq. 
Smaller nanostructures feature larger hybridization and 
are less prone to develop magnetic order. In the case of 
hexagons we expect a critical size above which exchange 
interactions take over and the edges magnetize. This is 
indeed what we have obtained from our mean field cal- 
culations. 

The local magnetization to/ for the N = 8 hexagon 
with U = 2.5eV is shown in the right panel of figure 3. 
The local magnetic moments lie mainly on the edges. We 
quantify the formation of local moments in compensated 
structures by the sublattice resolved average magnetic 
moment on the edge atoms 

{m^)edge = ^ (5) 

where r] = A, B and the sum runs over the 3iV edge atoms 
of the J] sublattice in the hexagon. For a given value of 
U, there is a critical value of N below which this quan- 
tity is zero. In Fig. [3^ we plot {niri)edge as a function of 
N for three different values of U. We always find that 
{m a) edge = ^{'niB)edge- This panel also shows how the 
critical size depends on U. When sweeping J7 in a rather 
wide range {U = 1.5 to U — 3.5 eV) the largest possible 
paramagnetic hexagon goes from 7 to 4. We have also 
estimated the critical size with the help of DFT calcula- 
tions and found that the largest paramagnetic hexagon 
corresponds to iV = 8, which is consistent with the mean 
field Hubbard results for small U = 1.5eV. 

Final remarks and conclusions.- We have seen how the 
magnetic properties of graphene nanostructures are in- 
timately related to the sublattice imbalance Na — Nb 
in agreement with Lieb's theorem^l. This is related to 
previous work on vacancies in graphene^^. As a con- 
sequence of Lieb's theorem a single vacancy results in 
the formation of a local moment with 5 = 1/2 and the 
sign of the spin coupling between two single atom vacan- 
cies would depend on whether or not they belong to the 
same sub-lattice-^. The correlation between sublattice 
and sign of the exchange interaction is also seen in our re- 
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FIG. 3: (Color online), (a) Sublattice resolved average mag- 
netic moment (eq. ((5]l) as a function of A'^ for 3 values of U. 
(b) Magnetization density for U — 2.5eV and N = 8. Arrows 
are plotted vertically for the sake of clarity 



suits for triangular and hexagonal nanoislands: Moments 
in the same sublattice couple ferromagnetically whereas 
moments in different sublattice couple antiferromagnet- 
ically. Indirect exchange interaction in graphene follows 
the same rule^i. 

Nanomagnets show remanence and hysteresys because 
of magnetic anisotropy, which originates in the spin-orbit 
interaction, very small in graphene. Therefore, the direc- 
tion of the spontaneous magnetization, M, of graphene 
nano-islands will fluctutate in the absence of an applied 
magnetic field. At zero field, the detection of magnetism 
should rely on properties that depend on |A/|, the mod- 
ulus of the magnetization vector. An example of this is 
the quasiparticle density of states, as probed with single 



electron transport in systems with spin polarization and 
without magnetic anisotropy^S. The controlled addition 
of single electrons to other nanomagnetic structures, like 
magnetic semiconductor quantum dots^^, afford the elec- 
trical control of their magnetic properties. This deserves 
further attention in the case of magnetic graphene nano- 
islands. 

In conclusion, we have studied the emergence of mag- 
netism in graphene nano-islands with well-defined zigzag 
edges. Our DFT calculations suggest that the mag- 
netic structure of the ground state of graphene nanois- 
lands can be described with a simple Hubbard model. 
Ground states with finite spin S appear in structures in 
which the number of atoms of one of the sublattices is 
larger than the other, Na > Nb, like triangular islands. 
The single particle spectrum of these structures features 
Nz = Na — Nb states with strictly zero energy, local- 
ized in the A sublattice, which yield a magnetic ground 
state with finite magnetic moment S = when 
interactions are included, both in a mean field Hubbard 
model and with DFT calculations. Compensated struc- 
tures {Na = Nb) like hexagons have 5* = 0. However, 
they develop spontaneous sublattice magnetization above 
a critical size. All our results are nicely consistent with 
Lieb's theorem^i and complement the theorem in the case 
of S" = ground states. 

Note added: Upon the completion of this work we have 
been aware of a related work by E. Ezawa-^^, in the single- 
particle approximation, and De-en Jiang et al. and O. 
Hod et al. doing DFT calculations^. 
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